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ABSTRACT 



Context. Little is known about the physical properties of cometary nuclei. Apart from space mission targets, measuring the thermal 
emission of a nucleus is one of the few means to derive its size, independently of its albedo, and to constrain some of its thermal 
properties. This emission is difficult to detect from Earth but space telescopes (Infrared Space Observatory, Spitzer Space Telescope, 
Herschel Space Observatory) allow reliable measurements in the infrared and the sub-millimetre domains. 

Aims. We aim at better characterizing the thermal properties of the nucleus of comet 8P/Tuttle using multi-wavelentgh space- and 
ground-based observations, in the visible, infrared, and millimetre range. 

Methods. We used the Plateau de Bure Interferometer to measure the millimetre thermal emission of comet 8P/Tuttle at 240 GHz 
(1.25 mm) and analysed the observations with the shape model derived from Hubble Space Telescope observations and the nucleus 
size derived from Spitzer Space Telescope observations. 

Results. We report on the first detection of the millimetre thermal emission of a cometary nucleus since comet C/1995 01 Hale-Bopp 
in 1997. Using the two contact spheres shape model derived from Hubble Space Telescope observations, we constrained the thermal 
properties of the nucleus. Our millimetre observations are best match with: i) a thermal inertia iower than ~10 J K"' nr^ s"''^, ii) an 
emissivity lower than 0.8, indicating a non-negligible contribution of the colder sub-surface layers to the outcoming millimetre flux. 

Key words. Comet: individual: 8P/Tuttle - Radio continuum: planetary systems - Techniques: interferometric 



1. Introduction 

Much of the scientific interest in comets stems from their po- 
tential role in elucidating the processes responsible for the for- 
mation and evolution of the Solar System. They appeared in the 

' outer regions of the protoplanetary disk, when the giant plan- 
ets were formed by core-accretion ('Safronov & Zviagina|[ 
Pollack etal. 1996) or disk instability (Cameron 1978t 

' [2003 ). Together with the asteroids. Centaurs, and transneptu- 

. nian objects, comets are the remnants of the planetesimals not 
accumulated into the planets. Comets formed in the giant plan- 
ets region were ejected to the outer Solar System and form the 
Oort Cloud. Some of them return to the inner Solar System as 

' long-period c omets (now called nearly isotropic comets NIC, 
iDuncanet a n il988.) . such as C/1995 Ol Hale-Bopp, or become 
periodic comets after several perihelion passages and shorten- 
ing of their orbit due to gravitational perturbations. The later 
ones constitute the Halley-family comets (or returning NIC) to 
which IP/Halley and 109P/Swift-Tuttle are significant represen- 
tatives. Comet 8P/Tuttle also belongs to this group: its orbit has a 
high ecliptical inclination (55°) with a short period (13.5 years). 
On the other hand, it is believed that the Kuiper Belt, gravi- 
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tationally influenced by the outer planets, supplied the popu- 
lation of Jupiter-family short-period comets (now called eclip- 
tic comets EC, Duncan et al. 1988). Comets 19P/Borrelly and 
9P/Tempel 1, investigated by the Deep Space 1 and Deep Impact 
missions respectively, are typical comets of this category. 

Given their different origins and the different evolution 
schemes followed since their formation, one could expect to ob- 
serve different physical properties for the NICs and ECs. For 
the time being, no obvious correlation between the chemical and 
physical properties of a comet and its dynamical class has been 
measured (Crovisier et al. 2009). 

One of the possible ways to investigate differences between 
the two classes of comets is to measure their size distribution. 
Reliable siz e determin ations have been obtained for only 13 
NICs (iLamv et al.1 '2004). using ground- and space-based tele- 
scopes in the visible and infrared domain. The range of radii is 
surprisingly broad, from 0.4 to 37 km, much broader than that 
of the ECs but with only 13 objects, robust conclusions can- 
not be drawn. Further measurements are required to confirm this 
trend. In this context, it is important to develop new methods 
to obtain reliable size estimates from the Earth. The millimetre 
wavelength range is well suited for this purpose because the at- 
mospheric transmission is better than in the infrared, and also 
because the ALMA observatory will offer unique observing ca- 
pabilities in the near future. However, the flux of a comet nucleus 
is much lower in the millimetre domain than in the infrared, and 
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observations are still challenging and reserved to a few bright 
comets. When available, data sets at different wavelengths (vis- 
ible, infrared, millimetre) complete one another and allow de- 
tailed studies of the nucleus properties. 

Comet 8P/Tuttle was observed in 1992 using the Nordic 
Optical Telescope (NOT) at 6.3 AU from the Sun and appeared 
inactive at this time. From its apparent magnitude, the nucleus 
diameter was estimated to 15.6 km (Licandro et al. 2000), but 
photographic observations performed in 1980 at 2.3 AU sug- 
gested a nucleus 3 times smaller or highly elongated. The close 
approach to the Earth in December 2007-January 2008 (0.25 
AU) and its supposed large size made it a very interesting target 
for Earth based observations. In 2007-2008, the comet was ob- 
served with different techniques in the visible (ligh tcurve mea- 
surements with the Hubble Space Telescope (HST). lLamv et all 
[2010b), the infrared (Spitzer Sp ace Telescope (SST) observa- 
tions, |GroussineLal. 2008,'2010 |) and th e radio domains (radar 
experiments with Arecibo. Harmon et alJl2010,) . To complement 
this multi-wavelength set of observations, we observed comet 
8P/Tuttle with the Plateau de Bure Interferometer to measure 
the millimetre thermal emission of its nucleus. The only similar 
measurement of a comet nucleus was pe rformed more than 1 
years ago, in 1997 for comet Hale-Bopp (lAltenhoff et alJll999h . 
which demonstrates the difficutly of such observations and the 
unique opportunity offered by comet 8P/Tuttle. 

In this paper we present the results of our observations of 
comet 8P/Tuttle carried out with the Plateau de Bure interferom- 
eter at 1.25 mm (240 GHz) in December 2007-Januai-y 2008. A 
description of the observations is given in Sect. IH W e analysed 
the data using a shape model of Lamv et alj ( 1201 Obi) and a ther- 
mal model of the nucleus, which are described in Sect. [3] The 
results are presented in Sect. |4] and summarized in Sect.|5] 

2. Observations 

Comet 8P/Tuttle has been observed twice at a wavelength of 
1.25 mm (240 GHz) with the IRAM interferometer on 29 
December 2007 and 8 January 2008 (Table©- The IRAM inter- 
ferometer is a 6 antennas (15 m each) array located at the Plateau 
de Bure, in the French Alps, and equipped with heterodyne, dual 
polarization, receivers operating around 1, 2 and 3 mm (230, 
150, and 100 GHz, respectively). On both observing dates the 
array was set in a compact configuration with baseline lengths 
ranging from ~20 to 150 m. This results in a synthesized beam 
diameter of about 1 .2" at 240 GHz. The whole calibration pro- 
cess was performed using the GILDAS software packages de- 
veloped by IRAM tPetVi,2005l) . 

29 Dec. 2007: For this first attempt, the comet was tracked 
using an ephemeris computed using the JPL solution K074/19 
for its orbital elements. The instrument was tuned in a way to 
search for CH3OH lines around 241 GHz: 4 narrow correlator 
units were dedicated to the CH3OH lines, and a total bandwidth 
of 1 .2 GHz was used to measure the continuum emission of the 
comet. 8P/Tuttle was observed from 15 h to 23 h UT, under cor- 
rect conditions (system temperature of 200 K and a phase rms 
of -30°) during the first half, until 19 h UT. The second half of 
the observations suffered bad conditions (unstable phase, cloudy 
weather) and is not usable. The final dataset represents ~1.6 h 
on source aquired between 15 h UT and 19 h UT. The gain 
calibration sources (observed every 22 minutes to monitor the 
instrumental and atmospheric phase and amplitude variations) 
were 0133-t-476 and 0059+581. MWC349 was used to deter- 



mine the absolute flux scale. The pointing and focus corrections 
were measured every 40 and 80 minutes, respectively. Some sin- 
gle dish spectra of the CH3OH lines were recorded all along the 
observing period. The sky contribution was cancelled in posi- 
tion switching mode (ON-OFF), with a OFF position distant of 
5' from the comet. We do not retrieve any CH3OH line detection 
in interferometric mode but the lines were detected in the ON- 
OFF spectra. Their analysis is presented in Biver et al. (2008i) . 
together with other single dish observations of 8P/Tuttle. 

Given the distance of the comet, the apparent diameter of 
a nucleus with a radius of 10 km is about 0.1" and we do 
not expect to resolve the nucleus of 8P/Tuttle in the Plateau 
de Bure data. As a result, we analysed the continuum observa- 
tions in the Fourier plane, fitting the Fourier Transform (FT) of 
a point source to the observed visibilities. We measure a flux 
F = 2.4 + 0.7 mJy located at an offset of (-0.5",-1.4") in (RA, 
Dec) with respect to the pointed position (the astrometric pre- 
cision, estimated dividing the beam size by the signal to noise 
ratio, is ~ 0.3"). Compared to the latest ephemeris solution (JPL 
K074/27, including observations performed up to October 2008) 
the offset is (0.4",3.8"). To iflustrate the result of this fit, we 
present in Fig. [1] the real part of the visibilities as a function of 
the uv-radius. The visibility table has been shifted to the position 
where the point source was found in the fit, so that the imaginary 
part of the visibilities is null and their real part is the source flux. 
On the figure, the real part of the visibilities has been averaged 
in 30 m bins in uv-radius. We overplot the point source flux that 
was found by the fitting procedure. The absolute flux calibration 
adds an uncertainty of 20% on both the flux and its uncertainty. 
We estimate the overall ±lcr limits as - 1.2 x{F + Icr) = 
3.7 mJy and F^''^ = 0.8 x(F- Icr) = 1.4 mJy. From this we de- 
duce the final flux and corresponding eiTorbars F = 2.4;^[ y mJy. 



08 Jan. 2008: An updated solution of orbital elements (JPL 
K074/21) was used to compute the ephemeris of the comet at 
this date. Our single dish observation s of the comet at th e IRAM- 
30 m telescope (in early January, iBiver et al ] l2008h showed 
that its gas production was not sufficient to enable interfero- 
metric study of the coma. As a result we dedicated all the cor- 
relator units of the Plateau de Bure to continuum observations 
for this second run, without changing the observing frequency 
240 GHz. This results in a bandwidth of ~2 GHz, which repre- 
sents a gain of a factor 1 .3 in point source sensitivity with respect 
to the December observations. The comet was observed from 
15 h UT to 21 h UT, under poor, degrading with time, phase 
stability conditions (a system temperature of 200 K and a phase 
rms of ~30-100°, depending on the baseline length). As a re- 
sult we use in the analysis only the 2 first hours of observations 
(1.3 h on source), when the conditions were the best. Although 
the integration time is longer for the December observing run, 
the larger bandwidth in January results in a lower thermal noise 
for this dataset. The gain calibration sources were 0235+164 and 
0048-097 and 3C454.3 was used as a reference to determine the 
absolute flux scale. 

Fitting the FT of a point source to the observed visibilities in 
the Fourier plane, we retrieve a flux of 3.0+0.5 mJy located at an 
offset of (1.4",-0.1") in (RA, Dec) with respect to the pointed 
position (the astrometric precision is ~ 0.2"). Compared to the 
latest solution (JPL K074/27) the off'set is (1.8", 1.2"). This re- 
sult is illustrated with the real part of the visibilities presented as 
a fuction of the uv-radius in Fig.[T] Taking into account the flux 
calibration uncertainty, we obtain a flux of 3.0;^[ q mJy (Icr) or 
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Date 


UT 
h 


Ephemeris 


A 

AU 


ri, 

AU 




Flux* 
mjy 




29 Dec. 2007 
08 Jan. 2008 


15-19 
15-17 


JPL K074/19 
JPL K074/21 


0.26 
0.28 


1.11 
1.06 


54° 
65° 


2.4+0.7 
3.0+0.5 



" Phase angle of the observations 

* The flux and its uncertainty are measured by a fit of a point source to the observed visibilities. The absolute flux calibration adds an uncertainty 
of 20%. 
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Fig. 1. Real part of the visibilities as a function of uv-radius for 
the December (top) and January data (bottom). The visibility ta- 
bles have been shifted to the position found by the fitting pro- 
cedure so that their real part represents the point source flux. To 
increase the signal to noise ratio on individual points, the visibil- 
ities have been averaged in 30 m bins of uv-radius. The result of 
the fitting procedure in the uv-plane is presented by the dashed, 
grey line. 



3.0;^j g mJy (3cr). In the following analysis, we only use the flux 
measured on January 8, as it offers a better signal to noise ratio. 

We assume that the detected emission is solely due to the nu- 
cleus thermal emission. The contribution of dust thermal emis- 
sion is expected to be weak, based on previous observations of 
the d ust continuum in the millimetre an d submillimetre domains 
(e.g. iJewitt&Luul 119921; iJewitt & Matthe ws 1997). For exem- 
ple, using the JCMT (HPBW = 19 "^ iJewit t & Matthews (Il997h 
measured a flux of 6.4 mJy at 350 ;um for comet C/1996 B2 
(Hyakutake) at r/, = 1.08 AU and A = 0.12 AU. Assumin g that 
the d ust opacity varies according to /I " dJewitt & Matthewsl 
Il997i) and that the dust brightness distribution varies according 
to p"', where p is the distance to nucleus, we derive a dust flux 
at 240 GHz of 0.05 mJy for a beam size of l"and a geocentric 
distance of 0.26 AU. The gaseous activity of 8P/Tuttle in early 
January 2008 was typically 5 -10 times lower tha t the activity 
of Hyakutake at r/, - 1.08 AU dMumma et al.lll996tlBonev et al.) 



l2008llKobavashi et al.l20Tol) . Hence, the contribution of the dust 
emission to the detected flux at 240 GHz in comet 8P/Tuttle 
should not exceed 0.01 mJy, compared to the 3.0 mJy of the 
nucleus. 

3. Nucleus thermal Model 

3.1. Shape model 



iHarmon et al.l (l2008l l2010h interpreted their Arecibo radar ob- 
servations of the nucleus of comet 8P/Tuttle as implying a 
contact binary and p roposed a shape mo del composed of two 
spheroids in contact. ' Lamv et aP (1201 Obi) found that the light 
curve of the nucleus derived from their Hubble Space Telescope 
observations indeed was best explained by a binary configura- 
tion and derived a model composed of two spheres in contact. 
However these two models profoundly differ first in the direc- 
tion of the rotational axis and second, in the size ratio of the 
two components of the nucleus. FGroussin et al.l (2010) used both 
shape models to interpret Spitzer Space Telescope thermal light 
curves of th e nucelus and concluded that the shape model of 
iLamv et al] (|2010b) bett er matches the SST observations. 

The shape model of Lamv et al. (2010b) consists of two con- 
tact spheres with respective radii of 2.6±0. 1 km and 1 . 1 +0. 1 km; 
a sphere with a radius of 2.8 km would have the same cross- 
section. The pole orientation is RA - 285° and Dec - +20°, 
which gives an aspect angle (defined as the angle between the 
spin vector and the comet-Earth vector) of 82° on 29 Dec. 2007 
and 103° on 08 Jan. 2008, close to an equatorial view. The 
rotational per i od is 1 1 .4 h dLamv et al.l 12008) in this model. 
iHarmon et al.l (12010 ') derived a more precise value of the pe- 
riod based on the radar experiments carried out at Arecibo: 
11.385+0.004 h. By rotating the nucleus back from the radar 
epoch to the HST epoch, we were able to connect the rotation 
phases and thus determine the true or siderea l rotation period 
P,id = 1 1.444+0.001 h (see 'Lamv et al.l l2010bl) . 

In view of the above discussion we used two different shape 
models for our analysis: i) a simple spherical shape model with 
a radius of 2.8 km and ii) the more complex and realistic shape 
model of iLamv et alJ (l20I0bl) . made of two contact spheres. In 
the later model, the shape is divided into 2560 triangular facets 
of comparable size. We consider the size for these two shape 
models as robust, so that it is a fixed parameter. 

3.2. Thermal model 

The interpretation of the millimetric observations requires a ther- 
mal model for the nucleus. We used the thermal model presented 
in Groussin et al. (2010) for comet 8P/Tuttle and alre ady ex- 
tensively described in s everal past articles (e.g., Gro ussin et al.l 
l2004l ; lLamv et alj|2010a ). For each facet of the shape model, we 
solve for the surface energy balance between the flux received 
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from the Sun, the re-radiated flux, and the heat conduction into 
the nucleus. As the nucleus rotates around its spin axis, the illu- 
mination changes, and the heat conduction equation is computed 
for each facet considering a one-dimensional time-dependent 
equation. The projected shadows are taken into account in our 
model. As a result, we obtained the temperature of each facet as 
a function of time, over one rotation period. We then integrate 
the flux over each facet to calculate the total millimetric flux 
received by the observer as a function of time and wavelength 
(1.25 mm). 

T he active area o f 8P/Tuttle is restricted to < 15 % of the sur- 
face (iGroussin et al. 2010). As for 9P/Tempel 1 with an active 
area of 9% (L isse et al., 2005) , the sublimation of water ice can 
be neglected in the energy balance for the calculation of the ther - 
mal flux emitted from the nucleus surface (IGroussin et al.l2007h . 



4. Results 

Using the above simple spherical shape model with a radius 
of 2.8 km and our thermal model with the parameters of 
IGroussin etaP (l2008l) (thermal inertia / = J K ' m ^ s '^^ 
beaming factor rj = 0.7, emissivity e - 0.95), which 
is in that case identical to the Standard Thermal Model 
(iLebofskv & Spenceil [T989I) . we obtain a millimetric flux of 
5.6 mJy at 1.25 mm on 8 Jan. 2008. This is larger than, though 
barely in agreement with, our observed flux of 3.0^[ g mJy (Icr) 
or 3.0^j g mJy (3cr), which in fact corresponds to a smaller radius 
of r = 2.0 + 0.4 km (Icr) using the same thermal model and pa- 
rameters. In order to investigate the origin of this discrepancy we 
studied the influence of several parameters in the calculation of 
the flux at 1 .25 mm: i) shape, ii) thermal inertia /, iii) beaming 
factor T], and iv) emissivity e from the surface and sub-surface 
layers. Our goal is to investigate which of these parameters can 
help to decrease the nucleus flux to make it compatible with our 
observations. 



4.1. Shape effect 

The simple calculation described above has been made using 
our spherical shape model with a radius of 2.8 km. As de- 
scribed in Sect. [3] a more reali stic shape model exists , com- 
posed of two contact spheres (Lamv et"ani2008l |2010b). With 
this shape model, the flux is no more constant with time as for 
the spherical case, but changes during the nucleus rotation to 
produce a thermal lightcurve. Under our standard assumptions 
(/ = 0, 77 = 0.7, e - 0.95) we calculated the synthetic thermal 
light curve of this shape model, as illustrated in Fig.|2| The flux 
varies with time, with a minimum of ~4.45 mJy and a maxi- 
mum of ~5.65 mJy. The maximum flux is reached when the pri- 
mary and secondary are both illuminated, and this is equivalent 
to a spherical nucleus of radius 2.8 km as explained in Sect. 13.1! 
The minimum flux is reached when the secondary eclipses the 
primary, and this is equivalent to a spherical nucleus of radius 
2.5 km. As the amplitude of the ligth curve is large (1.2 mJy), 
depending on the time of observation, the observed flux can be 
quite different. 

Using the sidereal rotation period of 11.444+0.001 h, we 
phased the observations performed at Plateau de Bure on 8 
Jan. 2 008 with the Arecibo radar experiments of iHarmon et al.l 
(I2OI OI) performed in early January 2008. The reference epoch 
of the Arecibo radar study is Jan. 4.0046, i.e. ~10 rotation peri- 
ods earlier than our observations, which translates to an uncer- 
tainty of 0.0004 day, i.e. less than 1 minute. At the reference 
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Fig. 2. Thermal lightcurve of 8P/Tuttle as seen from the Plateau 
de Bure on 8 January 2008. We used our standard set of pa- 
rameters for the thermal inertia (/ - 0), the beaming factor 
(rj = 0.7) and the emissivity (e = 0.95). The shape model is 
that o f 8P/Tuttle derived from H ubble Space Telescope observa- 
tions jLamv et al.ll2008ll2010bl) . The lightcurve has been p hased 
relatively to the radar observations of^ Harmon et alJ (l2010l) . per- 
formed ~10 rotation periods earlier. The maximum (resp. min- 
imum) flux corresponds to the emission of a spherical nucleus 
of 2.8 km (resp. 2.5 km). The horizontal thick solid line corre- 
sponds to the observation time at Plateau de Bure (2 h long). 



epoch, the rotation phase is -5+ 1 .5° with 0° being the time when 
the nucleus is seen b roadside, with the larger lobe approaching 
dHarmon et al.ll20l"ol) and corresponds to the maximum that fol- 
lows the lightcurve minimum in Fig.|2l on Jan. 8.75. The 1.5° 
error bar on the reference adds another 3 min of uncertainty on 
the rephasing. The total uncertainty does not exceed 5 min and 
is small in comparison with the length of the Plateau de Bure 
observations (2 h). 

The lightcurve presented in Fig.|2]has been rephased and one 
can see that the Plateau de Bure observations (around Jan. 8.67) 
were performed around its minimum. Under the standard param- 
eters, we expect a flux of ~4.45 mJy at that time. This is still 
higher than the 3.0;^[ q mJy measured with the Plateau de Bure, 
thus we assessed the impact of different thermal parameters of 
the nucleus to investigate this discrepancy. 

4.2. Thermal inertia effect 

We present in Fig. [3] the thermal lightcurve expected from the 
comet 8P/Tuttle for different values of the nucleus thermal in- 
ertia. When the thermal inertia increases, the diurnal temper- 
ature variations decrease: the temperatures become cooler on 
the day side and warmer on the night side. In our case, since 
the phase angle is large (65°), the heating on the night side is 
more pronounced than the cooling on the day side, and overall, 
the observed flux becomes larger as thermal inertia increases. 
It is clear from Fig. [3] that a larger thermal inertia does not 
help to solve our issue as the flux only gets larger. In the range 
0-200 J K~' rn~^ s~ '^^ for the thermal inertia estimated by 
iGroussin et al.l (l2008l) . we then favour the lower values (typically 
<10 JK-i m-2 
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Fig. 3. Expected thermal lightcurves of the comet 8P/Tuttle at 
1.25 mm for different values of the nucleus thermal inertia be- 
tween (in red) and 300 MKS (in blue). The beaming factor 
is 77 = 0.7 and the emissivity e - 0.95. The geometrical con- 
ditions are that of the observations of 8P/Tuttle on 08 January 
2008 (ru = 1.07 AU, A = 0.28 AU, phase angle = 65°). 
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Fig. 5. Expected thermal lightcurves of the comet 8P/Tuttle at 
1 .25 mm for different values of the surface emissivity e rang- 
ing from 1 (highest, solid red line) to 0.4 (lowest, dashed yel- 
low line). The thermal inertia is / = and the beaming factor 
77 = 0.7. The geometrical conditions are identical to that of the 
observations of 8P/Tuttle in January 2008. 
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Fig. 4. Expected thermal lightcurves of the comet 8P/Tuttle at 
1.25 mm for different values of the surface beaming factor (77). 
The thermal inertia is / = and the emissivity e - 0.95. The 
geometrical conditions are identical to that of the observations 
of 8P/Tuttle in January 2008. 



4.3. Beaming factor effect 

We present in Fig. |4]the thermal lightcurve emitted by our syn- 
thesized, bilobate nucleus as a function of the beaming factor 
T]. In our thermal model, the beam ing factor follows the strict 
definition given b v 'LagerrosI (1 1 998h where it only reflects the in- 
fluence of surface roughne ss and ?? must b e lower or equal to 1 .0. 
In addition, according to lLagerrosI (Il998h . rj must be larger than 
0.7 to avoid unrealistic roughness, with r.m.s. slopes exceeding 
45°. The acceptable range for rj is then 0.7-1.0, as illustrated in 
Fig.|4] When roughness increases (i.e., t] decreases), the surface 
temperature increases due to self-heating on the surface, and the 
total flux increases. Changing r] from 0.7 dGroussin et alj|2008h 
to 1.0 decreases the minimum flux from 4.45 mJy to 4.1 mJy, 



which is a small improvement but not sufficient to explain the 

observed flux of 3.0 mJy. 

T he SST observations of comet 8P/Tuttle dGroussin et al.l 

I2008ll2010t) suggested a beaming factor rj on the order of 0.7. 
This value is constrained by the infrared spectrograph (IRS) ob- 
servations, which lasted only 10 minutes, a short time compared 
to the rotation period of ~11.4 h. The SST observations per- 
formed on 2.76 Nov. 2007 and, 141 periods later, the nucleus 
was in the same rotation phase on 8.65+0.02 Jan, 2008, exactly 
at the time of the Plateau de Bure observations (Fig. |2]l. As a 
result, although higher t] values correspond to lower millimetre 
fluxes, we favour the rj - 0.7, which is well constrained by the 
Spi tzer ob servations. 

lDeib6etal.l (120031) suggested an increase of the beaming 
factor with phase angle for Near-Earth Asteroids (NEAs). If 
confirmed, this trends could explain why our observations per- 
formed at 65° phase angle favor a larger value of 77, close to one, 
while SST observations performed at 39° phase a ngle favor a 
smalle r value, close to 0.7. However, according to Delbo et al.l 
(l200l . their statistic is low and more observations are required 
to confirm this trend. 

Finally, it cannot be excluded that the effect of the rough- 
ness (77) may be wavelength dependent. Roughness at micron 
scale could indeed be more important than at millimetre scale, 
resulting in a smaller beaming factor in the infrared. But this 
point cannot be confirmed since we cannot rule out a fractal sur- 
face, in which case the roughness would be more or less scale- 
independent. 

4.4. Emissivity and depth effect 

The thermal flux emitted from the nucleus is proportional to 
the surface emissivity, as illustrated in Fig. |5] The flux at the 
lightcurve minimum is in Icr agreement with our observations 
for an emissivity lower than 0.8. While the surface emissivity 
of small bodies in the mid-infrared is close to unity, it is still 
unknown at millimetre wavelengths and might be lower than 
in the mid-infrared. This point has already been addressed by 
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iFernanded (l2002h for comet C/1995 Ol Hale-Bopp, where an 
emissivity of 0.5 is required in the milHmetic range to reconcile 
the nucleus sizes derived from infrared and millimetric obser- 
vations. Based on a steroid observations in the sub-millimetre, 
iRedman et alj (Il998 h found that the emissivity of 7 asteroids de- 
creases with increasing wavelength. The observations of asteroid 
(2867) Steins at 1.6 and 0.53 mm performed with the radiome- 
ter MIRO o riboard the Rosetta spacecraft confirm this trend: 
iGuUds et al.l ( 12010 ) found an emissivity decreasing from 0.85-9 
at 1.6 to 0.6-0.7 at 0.53 mm. A low milli metric emissivity of 
0.6 wa s also found for asteroid 4 Vesta by iMueller & LagerrosI 
([T998h . Overall, this means that our values (e < 0.8) are plausi- 
ble. 

A possible physical explanation for the lower emissivity is 
that a part of the millimetre thermal flux arrises from sub-surface 
layers which are colder than the surface itself. Such a tempera- 
ture gradient is expected given the low thermal intertia of the 
nucleus. The depth of the main contributing layer depends on 
the material opacity, which is unknown. Measurements on rocks 
showed that in most materials this depth lies between 3 and 100/i 
dCampbell & Ulrichs 1969). A depth of lOA is commonly used 
when dealing with plan etary surfaces (e.g. 11/t is used for the 
martian surface regolith. FMuhleman & Ber^ell991l;lGoldin et al.l 
Il997h . Depending on the upper layers temperature gradient (con- 
troled by the thermal inertia) and the relative contribution of 
the different layers to the overall emission, the resulting spec- 
tral emissivity of the surface can be significantly lower than the 
usual infrared value close to unity. 

In order to investigate this effect, we calculated the nucleus 
integrated flux as a function of depth (down to 20 cm), relative 
to the surface flux, for different thermal inertia values. Figure |6] 
illustrates our results: as we observed the afternoon side of the 
nucleus, the shadowed part that we see is still warm because 
of thermal inertia and the contribution of the first sub-surface 
layers is comparable or larger to the surface, down to a certain 
depth (that depends on the thermal inertia) where it drops since 
diurnal temperature variations become negligible. 

Assuming a Poissonian weighted function (that peaks at 10 A) 
to describe the relative contributions of the sub-surface layers, 
w hich seem s to be a good approximation compared to the work 
of lMoulletl(l2008ft . we integrate the weigthed flux from Fig.|6]as 
a function of depth, to derive the effective outcoming flux com- 
pared to a pure surface emission. The results are presented in 
Table |2] Depending on the thermal inertia, the flux is reduced 
by 9% to 17% compared to a pure surface emission. For low 
thermal inertia (<10 J K"' m"^ s"'''^), which are in better agree- 
ment with our observations as explained in Sect. 14.21 the flux is 
reduced by 14-17%, which coiTesponds to an "effective" emis- 
sivity of 0.79-0.82 assuming our standard surface emissivity of 
0.95. According to Fig.|5]and as explained above, this allows to 
match our observations within Icr. 

Our results confirm that a low emissivity in the millimetric 
wavelength range can result from the emission of sub-surface 
layers with a low thermal inertia. In the mid-infrared, typically 
around 10 fim, 10 A corresponds to 0.1 mm, a depth over which 
the flux is close the surface flux (Fig. |6]l, explaining why the 
above effect only matters in the millimetric wavelength range 
and beyond, in the centimetric and radar wavelength range. 
However, for the centimetric and radar, there exists a few no- 
table exceptions, and in particular Ganymede, which centrimet- 
ric emission coiTesponds to a brightness temperature of about 
55-88 K, below a ny acceptable sub-surface temperature for this 
body dMuhleman & Berge 1991 ), and whose r adar reflection is 
highly anomalous (.Ostro & Shoemaker! 1 1990h . In this case, a 



Table 2. Nucleus thermal flux integrated over depth, compared 
to pure surface emission (1.0). 



Thermal inertia Relative depth-integrated flux 



3 0.83 

5 0.84 

10 0.86 

50 0.88 

100 0.91 



1.10 r ^ 

: 10 A 




: 1=100 : 
0.75 t : 

0.01 0.10 1.00 10.00 100.00 

Depth (mm) 

Fig. 6. Thermal emission from a sub-surface layer (relative to 
surface emission) as a function of its depth different values of 
the thermal inertia. 
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Fig. 7. Same Figure as Fig. |2]but for / = 0, e = 0.8, tj = 0.7 
(red curve) ox rf - 1.0 (blue curve). For rj - 1.0, which is our 
prefered value, about half of the thermal lightcurve agrees with 
our observations at the Icr level. 



higher reflectivity must also be invoked, likely resulting from 
backscattering of the surface. 
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5. Summary 

We report here the first detection of the thermal emis- 
sion of a cometary nucleus at millimetre wavelength since 
C/1995 Ol Hale-Bopp, performed at Plateau d e Bure on 8 Jan. 
2008. Using the rotation period derived bv lHarmon et alj(l2010i) . 
we phased our observations with respect to their radar exper- 
iments and found that Plateau de Bure observations were per- 
formed at a minimum of the lightcurve. We used the shape model 
developped by Lamv et al. ( 20 1 Ob) for 8P/Tuttle nucleus to com- 
pute thermal emission lightcurves using different thermal param- 
eter sets. From our analysis, we can conclude that to match our 
observed flux of 3.0;^} ^ mJy at 1.25 mm, the thermal inertia 
of the nucleus should be low, typically <10 J K"' m"^ s"'^^, 
in agreement with the ran ge 0-200 J K~ ' nT^ s ''^ derived 
from Spitzer observations dGroussin et alj 120081) . In addition, 
the millimetric emissivity of the nucleus should be lower than 
0.8. According to our thermal model, such a low value can be 
due to the non-negligible emission of the sub-surface layers of 
the nucleus, which are colder than the surface at depth of a 
few millimetre or more for low thermal inertia values. Similar 
and even lower values of the millimetric emissivity have been 
alrea dy mentionn ed to explain t he low fl ux emitted by aster- 
oids (^edman^etal. 1998; Muell er & Lage rros 1998; Fernandez 
I2OOI . As to the beaming factor, the lowest flux of the lightu- 
curve is in Icr agreement with our observed flux for any value 
in the plausible range 0.7-1.0, with a better agreement towards 
higher edge. However we favour rj = 0.7, the only value in agree - 
ment with Spitzer (IRS) observations dGroussin et alJ l2010l) . 
From the above conclusions, we generated a synthetic lightcurve 
for the nucleus of comet 8P/Tuttle, as observed from the Plateau 
de Bure on 8 Jan. 2008, using the following parameters: / = 0, 
6 - 0.8 and 77 in the range 0.7-1.0 to cover the possible values. 
The result is illustrated in Figure Q 

In the future, the Rosetta spacecraft, equipped wit h a mi- 
crow ave instrument operating at 190 and 562 GHz (Guil ds et alJ 
12007 ). will provide an unequaled dataset to study the surface 
of a cometary nucleus and its thermal properties. It will be also 
important to repeat ground-based studies on other comets, to de- 
termine whether all comets present similar properties or if they 
change from one comet to another, revealing different surface 
natures. The (sub-)millimetre interferometer ALMA will offer 
a significant gain in sensitivity and allow similar studies on a 
larger sample of comets, from both the ecliptic and the Oort 
cloud dynamical classes. 
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